Superfluid transitions in bosonic atom-molecule mixtures near Feshbach resonance 
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We study bosonic atoms near a Feshbach resonance, and predict that in addition to a standard 
normal and atomic superfluid phases, this system generically exhibits a distinct phase of matter: 
a molecular superfluid, where molecules are superfluid while atoms are not. We explore zero- and 
finite-temperature properties of the molecular superfluid (a bosonic, strong-coupling analog of a BCS 
superconductor), and study quantum and classical phase transitions between the normal, molecular 
superfluid and atomic superfluid states. 
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Experimental realizations and coherent manipulation 
of trapped degenerate gases 0, 0] is leading to excit- 
ing possibilities for studies of quantum liquids in previ- 
ously unexplored (e.g., extremely coherent and nonequi- 
librium) regimes. Magnetic field- induced Feshbach reso- 
nance (FBR) in ultracold atom collisions allows fine tun- 
ing of interactions in these quantum fluids, and was re- 
cently used to create a degenerate mixture of coherently- 
coupled alkali atoms and their diatomic molecules [3|. 
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FIG. I: Mean-field phase diagram for a bosonic atom- 
molecule mixture, showing molecular (MSF) and atomic 
(ASF) superfluid phases. 

In this Letter we study phases and phase transitions 
that take place in bosonic atom-molecule mixtures. Our 
main contribution is the prediction of a thermodynami- 
cally distinct "molecular superfluid" (MSF) phase, that, 
as illustrated in Figs. [21 ubiquitously intervenes be- 
tween the "normal" (N) and "atomic superfluid" (ASF) 
phases. Molecular superfluidity [and accompanying off- 
diagonal long-range molecular order (ODLRO)] distin- 
guishes MSF from the normal state, and the absence 
of atomic superfluidity from the ASF, in which both 
bosonic atoms and molecules display ODLRO. If atomic 



and molecular components can be imaged independently 
4], in a harmonic trap MSF should be easily identifi- 
able by a sharp Bose-Einstein condensation (BEC) peak 
in the molecular density profile and a broad, seemingly 
normal, thermal atomic cloud. 

As a conventional superfluid, MSF is characterized by 
a (molecular) acoustic second-sound mode. However, 
MSF also exhibits a gapped, Bogoliubov-like mode, de- 
rived from unpaired atom excitations. MSF ground state 
(bosonic analog of the BCS state) exhibits strong (atom 
and molecule) pairing correlations that in a trap should 
be observable in the atomic density-density correlation 
function. Experimentally, MSF should be accessible by 
tuning temperature, atomic density (or number), and de- 
tuning v. The MSF-ASF transition is in the (d + 1)- 
and d-dimensional Ising universality classes for T = 
and finite T, respectively, and is reentrant as a func- 
tion of detuning v and density n. The tricritical point, 
where N, MSF and ASF meet, exhibits nontrivial and, 
to our knowledge, unexplored quantum critical behavior 
for d < 4. We now sketch derivation of these results. 

Near a FBR a bosonic atom-molecule system is char- 
acterized by the grand-canonical Hamiltonian H fl = H — 
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(1) 



where (x) , ty a (x) are bosonic field operators for atoms 
(cr = I) and molecules [a = 2), h a = -(S 2 /2m ff )V 2 — fi a 
are the corresponding single particle Hamiltonians (fo- 
cusing for concreteness on the case of a homogeneous 
trap) with effective chemical potentials \i\ = [i and 
fi2 = 2/i — v. Chemical potential fi tunes the aver- 
age total number of atoms (whether free or bound into 
molecules) to N, and detuning v is related to the en- 
ergy of a molecule at rest, that can be experimentally 
controlled with a magnetic field. In the dilute gas limit 
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gi,g 2 ,gi 2 are proportional to the 2-body s-wave atom- 
atom, atom-molecule and molecule-molecule scattering 
lengths, respectively, and a characterizes coherent atom- 
molecule interconversion rate, encoding that molecules 
are composed of two atoms 0. 

The mean-field phase diagram as a function of 
^1,2 and f3 = can be worked out by min- 

imizing the imaginary-time (t) coherent-state action 

S = f rfr/AELlC^i + ff.ffti)]. Sim- 
pie analysis leads to three thermodynamically distinct 
phases (Fig. [TJ: (i) "normal" (N): *i = (^i) = 0, 
*20 = $2) = 0, (ii) "molecular superfluid" (MSF): 
*io = 0,1-20 ^ 0, (iii) "atomic superfluid" (ASF): 
"fio 7^ 0, ^>20 ^ 0. Condensed atoms cause a to act as an 
effective field on the molecular order parameter ^ 2 q, so 
an equilibrium phase in which atoms arc condensed, but 
molecules are not, is forbidden Q. 

We now examine in more detail these phases and cor- 
responding phase transitions. Phase N is stable for 
Mi, 2 < 0, with [i determined by the total atom constraint 
n = ni + 2ri2, which in the non-interacting limit, ap- 
propriate to a dilute weakly interacting gas, is given by: 



A d 



h/2 (e^)+2( rf + 2 )/ 2 / d/2 {e^- v) 



(2) 



where Ax = h/ ' \J2nm\kBT is the thermal de Broglie 
wavelength and f a {z) — Y^=i z n /n a {\ z \ < 1) is the 
extended zeta-function. 

The N-ASF transition line T c i(n, v) occurs at /i = for 
v > 0, while the N-MSF line T c2 (n, v) occurs at \i-2v = 
for v < (see Fig.[3J|. Using the appropriate asymptotics 
of fa(z), one obtains from (J2J: 
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baC 2 / d T c0 , 



\v\ < k B T c 



<•() 



\v\ > k B T c0 , 
(3) 

with c = i + 2( d + 2 )/ 2 , oi = 2( d+4 )/ 2 |r (2^) |/dcC(d/2), 

a 2 = 2- d ai , b x = 1, b 2 = 2-( d+2 )/ d and T c0 = 
(h 2 /2irm\kB)[nl 'c((d/2)] 2 / d the transition temperature 
at the tricritical point v = 0. 

In the neighborhood of T cl the "massive" molecular 
field ■tp 2 decouples at low energies (can be safely inte- 
grated out of the partition function, leading to an effec- 
tive quartic coupling g\ — > g~\ = gi — 2a 2 / 1 ) , and 
the N-ASF transition is identical to that of a single- 
component system, continuous so long as g\ > 0. At 
T = 0, the N-ASF transition takes place at vanish- 
ing atom density (the N phase is simply a vacuum of 
atoms), and although nontrivial, is exactly soluble 0, 
corresponding to a build-up of atomic superfluid (with 
condensate density nio = I'J'iol 2 ~ | A* 1 213 > with mean- field 
result (3=1/2 for d > 2, and j3 = d/4 for d < 2) as the 
trap is loaded. At T ^ the N-ASF transition lies in the 
usual d-dimensional XY-universality class 



Similarly, in the neighborhood of T c2 , ipi decouples and 
the resulting N-MSF transitions are in the same univer- 
sality classes discussed above. The full phase boundary is 
illustrated in Fig|21 In 3d it exhibits a square-root singu- 
larity at the tricritical point and for — > oo asymptotes 
to the single-component BEC temperatures T° 
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FIG. 2: Phase diagram for a bosonic atom-molecule mixture 
in d — 3, expressed in terms of detuning v and temperature 
T. It illustrates a finite-T tricritical point at T c o and a quan- 
tum critical point at u c (n, 0). In the weakly-interacting limit 
appropriate to experiments the ratio /T& = 2 5 ^ 3 . 



To study the MSF phase, we separate i/v = ^ a Q + 
<\> a into classical condensate fields "J CT o (with \&io = 
inside MSF) and fluctuations about it. Within MSF it 
is sufficient to expand to second order in fluctuations 
4> a , which leads to = £ (0) (* 2 o) + # (2) with: 
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The linear term in 



Ai = -2a*20, and A 2 = g 2 ^> 20 - 
4> 2 , (j>2 vanishes automatically by the self-consistent choice 
of ^20 as the true minimum of the free energy. To lowest 
order this gives 

n 2a = |*2o| 2 = (2n-u)/g 2 , (5) 

which coincides with the minimum of E^ and allows us 
to eliminate fi in favor of v and n 2 Q. 

For a homogeneous system, may be diagonal- 

ized by Fourier transformation 4> a = V~ x l 2 ^ k e ak x a CT k, 
followed by independent Bogoliubov transformations on 
atoms and molecules to new boson operators 7cr,k,7CTk : 
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k 7cr,k - |u<r,k! 



(0) 
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in which ev.k = £o-.k 



AV + r<r, £<T.k = h k 2 /2m„ and 



/ er k — |A CT | 2 , and iv,k, ^ have the same phase 
as X a . Using (J3J) one obtains 



\J [ei,k - f/2 - (.92/2 - gi2)n 2 o] 2 - 4a 2 n 2 o , 



^2,k = y £ 2,k + 2^2^20 £2,k ■ 



(7) 



The MSF ground state is defined as vacuum of atomic 
and molecular Bogoliubov quasi-particles 7 CTi k|MSF) = 0. 
It is easy to show that it is given by: 



I MSF) 
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10), 



(8) 



which is a coherent state in the k = molecular mode 
and shows BCS-like pairing correlations between time- 
reversed (k, — k) atomic and molecular single particle 
states. The amplitude Xo-,k = ?V,k/M<7,k is the Fourier 
transform of the wavefunction of k ^ atom (a = 1) 
and molecule (a — 2) pairs. 

As in a single-component SF the relation (JSJ (and more 
fundamentally, the Goldstone theorem |1C|1 ensures that 
MSF exhibits a gapless sound mode E 2 ,k ~ tiv k, for 
small k, with vq = yg^nw/2mi, corresponding to col- 
lective, long wavelength oscillations of the molecular con- 
densate. The resulting |k| singularity in ^2.k induces 
a long ranged power law tail in the Fourier transform 
X%( r ) ~ l/r d+1 . In addition, we find a gapped excitation 
branch describing atomic-like excitations, whose 

spectrum and eigenmodes (71 it, 7! k do not even carry 
a definite atom number) are qualitatively distinct from 
unpaired atomic excitations, ai^jO-fk of the N phase. 
The condition of being in the MSF is that the atomic 
gap £gip = y/e+e-, where 

e± = —v/2 - (52/2 - 312)^20 ± 2ay / nw, 



0, (9) 

remains positive. Correspondingly, the atomic pair wave- 
function xi( r ) ~ e~ r ^ 1 decays exponentially at large 
distance. The correlation length £1 = h/y / 2me- charac- 
terizes atom-pair size, and diverges as iSgap vanishes. The 
MSF- ASF transition takes place when the size of atom- 
pairs becomes comparable to the intermolecular separa- 
tion, and £1 is the diverging size of the coherent exchange 
loops of atoms between overlapping pairs. 

In the present quadratic approximation, the free en- 
ergy is given by 
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The condition of fixed density n = —(df/dfi) n20 , 



(10) 



(11) 



where the first term under the summation represents the 
T = interaction-induced condensate depletion, may be 
used to determine the condensate density n2o(n, v, T). At 
T = 0, in the weakly-interacting limit {o?n <C 1), in 3d 
we find: 
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d/2 



2 rf+ir(|) 



(12) 



where a\ 2 = a 2 — 2ai 2l a 2 ,ai 2 are molecule- molecule 
and atom-molecule scattering lengths, respectively, de- 
fined by (712,2 = 47rft 2 ai2.2/'7i2, the "detuning density" 
is n v — mi\v\/(Trh 2 ai 2 ), a = 2 3 / 2 mi\a\/(Trh 2 a~i 2 ), and 
X(y) = y 2 dxx d / 2 [(l + x) 2 — y 2 ]~ 3 / 2 is a scaling func- 
tion describing additional molecular-condensate deple- 
tion due to atom pairs. Equation (|12[) . together with the 
vanishing of e_, determines the MSF- ASF phase bound- 
ary. At T = 0, approximating n 2 Q ps n/2, one finds: 



(13) 



v c {n) w ~{g 2 /2 ~ gi 2 )n - 2u\J2n 



which is illustrated for 172 > 2gi2 in Fig. ^ At finite 
T, in the noninteracting limit, v a ^ — and u a k = 1, 
molecular condensate density reduces to 



T 

T'OC 
^9 



3/2 



2A?p 



/ 3/2 (e-l"l/»^) , 



(14) 

with the standard BEC result (first term) corrected by 
depletion due to thermal atomic excitations (last term), 
that is exponentially small at large \v\ and low T. 

From the structure of the interconversion (a) term in 
Hp it is clear that it is a Z 2 (Ising) symmetry that is 
broken at the MSF— + ASF transition [llj. We now show 
that for a homogeneous trap, both the T = and fi- 
nite T MSF- ASF transitions lie in the Ising universality 
class. This can be most easily seen from the coherent 
state action, which, when expressed in terms of real and 
imaginary parts of the atomic field ip\ = ipji + iipi, in 
MSF phase reduces to: 



S = 



dT / d d 



- 2iip I hd T ->p R + i/) R (hi - 2a^2o)ipi 



(15) 



where S- m t are terms not essential for our argument. In 
this form it is clear that in the presence of the molecular 
condensate, | 1 > 0, positive a reduces the 0{2) sym- 
metry to Z 2 , with %1>r reaching criticality before ipj. Be- 
cause canonically conjugate field tpi remains "massive" at 
the critical point [defined by where the coefficient of t\;\ 
vanishes, consistent with (|13fl ], it can be safely integrated 
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out and leads to a d + 1-dimensional (Lorentz-invariant) 
action even in the scalar order parameter ipR. Therefore, 
as asserted above, the T = MSF-ASF transition is in 
the (d-fl)-dimensional Ising universality class. The Ising 
transition is well studied, and leads to the following pre- 
dictions [T^j. F° r d = 3, up to logarithmic corrections, 
the mean-field theory derived above will be accurate. On 
the other hand in 2d, MSF-ASF exponents are nontriv- 
ial but are well-known. For example, standard scaling 
arguments predict: 

n lo ~W-v c \ 201 , EW~\v-v c \'"", (16) 

where /3j ss 0.31, zj = 1, and Vi ~ 0.63 are 3d Ising ex- 
ponents. These, together with the relevance of T at this 
quantum critical point also imply a universal shape of the 
MSF-ASF phase boundary u c (n,T) ~ v c (n,0) + aT 1 / u ' 
in Fig. [21 One may hope that when long-lived molecular 
condensates are produced, nontrivial behavior of Eg^ p (v) 
maybe observed in Ramsey fringes experiments 

Although from symmetry point of view ASF state is 
quite conventional, it also exhibits a set of interesting 
features not found in a single-component SF. For exam- 
ple, because it is a discrete (Ising) symmetry that distin- 
guishes ASF from MSF, it will exhibit a gapped (Ising) 
mode in addition to the gapless second-sound mode. As 
discussed above for a ^ 0, finite ^20 is always in- 
duced within the ASF. However, in the limit of small 
a, we predict a sharp crossover (indicated by a dashed 
^2 ff = — MiffWffl = curve in Figs. where 
the ip2 held would order for a = jfj) between small 
molecular condensate n2o ~ («Mi/5iA*2 ff ) 2 I0r ^ < 
and large n 2 o ~ /u| /fla f° r > 0, and vanishes as 
n-20 ~ (anio) 2 / s for — and weak a, with S a uni- 
versal XY-model magnetic field exponent (Sd=3 — 5) 0] . 

Another experimentally interesting feature of ASF is 
that a seemingly standard 2tt vortex in atomic conden- 
sate ^10 will generically split into two it vortices (see Fig. 
|3J) confined by a domain wall of length 

r Q = \J h 2 n% 2 /2manl , (17) 

that diverges as ASF^MSF phase boundary is ap- 
proached |l2j. This arises because in large a limit a 

71 7C 
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FIG. 3: 2ir atomic condensate vortex in ASF splits into n + ir 
vortex pair connected by a "normal" domain wall. 

2-7T vortex in atomic condensate induces a Air vortex in 



a molecular condensate. Such double molecular vortex is 
unstable to two fundamental 27r molecular vortices, that, 
in 2d repel logarithmically, but are confined linearly in- 
side the ASF phase. 

Finally, standard renormalization-group analysis 
shows that upper-critical dimensions for the tricritical 
point (dominated by scaling of a) are a%>° = 6 and 
d^=° = 4, and therefore should display nontrivial critical 
properties [l2j ]. 

We plan to present analysis of these problems as well as 
generalization to experimentally relevant harmonic trap 
in a future publication. 

To conclude, we have studied a mixture of bosonic 
atoms and their diatomic molecules and predicted an 
existence of a molecular supcrfluid phase, qualitatively 
distinct from the normal and atomic superfluid states. 
Ising transition MSF-ASF should be observable in such 
systems as a function of detuning u, temperature and/or 
density. 
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